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Abstract 

O ■ 

o : 

The equation of state for the pion gas is analyzed within the third virial 
approximation. The second virial coefficient is found from the 7T7r-scattering 
<]} \ data, while the third one is considered as a free parameter. The proposed 

Q | model leads to a first-order phase transition from the pion gas to a more dense 

t^j- ' phase at the temperature T p t < 136 MeV. Due to relatively low temperature 

this phase transition cannot be related to the deconfinement. This suggests 
. that a new phase of hadron matter — 'hot pion liquid' — may exist. 

>: 

^ | The widely used thermal models of the hadron gas (HG) are mostly based on the ideal 

gas approach. To take into account attractive interactions between hadrons, the hadron 
resonances are included in the ideal gas picture Q. This method, describing quite satis- 
factory the gaseous phase, does not lead to a phase transition (PT). In phenomenological 
models, PT is usually introduced 'by hand': it is postulated from the very beginning that 
two different phases, for instance HG and quark gluon plasma, exist; each phase has its own 
equation of state and Gibbs criterion is used to find the PT point. 

It is interesting to study whether a PT can naturally follow from the equation of state 
of real HG without any ad hoc assumptions. This would allow us to describe, at least at a 
>- ! qualitative level, both phases within a unique approach and to estimate the PT parameters 
^ ' using available hadron data. 

In the present letter we make the first step in this direction: we study the equation of 
state of the pion matter, i.e. we ignore (anti-)nucleons and strange particles. 

Our starting point is the equation of state for a gas of interacting particles in the form 
of virial expansion : 

oo 

p(T,n)=TY / a l (T)n l . (1) 
i=i 

It allows to calculate the gas pressure p at given temperature T and particle density n, 
provided that the virial coefficients «i(T) are known. 

In the case of a relativistic system, like a hot pion gas, the particle number is not con- 
served, hence it is inconvenient to choose n as an independent variable. On the other hand, 
the chemical potential \i for pions in isotopically symmetric systems at chemical equilibrium 
is equal to zero. Therefore, it is more convenient to choose T and [i as independent variables: 
n = n(T,fj, = 0) should be found from Eq.(ffl) in thermodynamically consistent way. 

Let us express the thermodynamical quantity of the gas satisfying Eq.(|I]) via T and [i. 
The particle density is given by the formula 
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Using the equation of state (0) leads to a differential equation for n(T, //): 

dp(T,n)dn(T,n) ^. ^ dn(T,n) 
n( - , ^ = ~d^ fy]— = T 2^ m i( T ) n ( T ^> 

Integration of this equation yields 
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where C(T) is the integration constant. In the low density limit n — > at fixed temperature 
(which is equivalent to — ► — oo) the equation (|j) should reproduce the properties of the 
ideal classical gas, i.e. n(T,fj) should coincide with the absolute activity z(T,fi), which is 
given by the formula Q: 
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where g is the (iso-)spin degeneracy factor (g = 3 for the pion), m is the pion mass, and 
4>(T; m) can be expressed via the K 2 modified Bessel function 
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This allows to fix the integration constant. Taking into account that ai(T) = 1 one gets 



C(T) 



exp(^; 



z(T,fi) g<t>(T;m) ' 
Finally, one has the following transcendental equation for n(T, /j,) 



(7) 



z(T, /i) = n(T, fi) exp 
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To apply this equation to the pion gas it is necessary to know the values of the virial 
coefficients. The second virial coefficient is straightly related to the second Mayer's cluster 
integral a 2 = ~ b 2 , which is given by 



& 2 = 4 0) +6«, 



where the first term 
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2 2g 2gm 2 T [K 2 {m/T)} 2 
appears due to Bose effects. It exists in an ideal Bose gas too. The second term 
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is related to the interactions between particles. The function Sl,i{b) stands for the phase 
shift of 7T7T-scattering in the state with angular momentum L, isospin / and center-of- 
mass energy e. The sum over L runs over those values, which ensure symmetry of the total 
two pion wave function, i.e. 

fO.2,4,6,... for 7 = 0,2 
I 1,3,5,7, ... for I = 1. 



The second cluster integral b 2 was calculated in |J. The dependence of a 2 on the temperature 
is shown in Fig.|TJ (see Ref. for details). 

A model independent determination of all other virial coefficients a,(T), i > 2 is hardly 
possible. The quantum theory is developed only for a 3 (T) flQfl , but lacking of detailed 
experimental information on the three-pion interaction makes it impossible to use this theory 
for the case of pion gas. For the higher virial coefficients, the problem looks even more 
hopeless. 

A simple model will be adopted for the further consideration: 

• The virial series is truncated after the third term: Oj(T) = 0, i > 3; 

• The third virial coefficient 03 (T), at fixed T, is considered as a free parameter. 

It is convenient to introduce the notation £(T) = 03 (T) / '[c^T)] 2 and treat the dimen- 
sionless quantity £(T) as a free parameter. In this case, the equation of state @ takes the 
form 

p(T, n(T, n)) = Tn(T, {l + a 2 (T)n(T, + £{T) [a 2 (T)n{T, fi)] 2 } , (13) 
where the value of n(T,fi) is found from the transcendental equation 

z(T, fi) = n(T, fi) exp 



2a 2 (T)n(T,fx) + ^(T) [a 2 (T)n{T 



(14) 



Let us analyze the behavior of the r.h.s. of Eq.flUP at some temperature T, as a function 
of n > 0. For convenience we denote it by f(n). Three distinct cases are possible (see Figs. 
0and|): 

1- £(T) — 0. The function f(n) has a single extremum — a maximum at the point 

i ^ *cn=o 

n a = <J (15) 

3|f(T)||aa| 11 J ^ u > 



2. ^(T) > 1/3. The function f(n) grows monotonically for all n > 0. 
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3. < £(T) < 1/3. In this case the function f(n) has two extrema: a maximum at the 
point 



1 ~ V 1 ~ 3g(T) 
n ° = J m i I ( 16 ) 



and a minimum at the point 



1 + .^1-3£(T) 



In the first case, the values of f(n) are bound from the above: f(n) < f(n a ). As z(T, fi) 
is a monotonically increasing function of the temperature, the equation ( |14|) has no solution 
at high temperatures T > T a , where the temperature T a is found from the equation 

z(T a ,fi) = f(n a ). (18) 

The numerical calculations give T a < 128 MeV assuming chemical equilibrium (/i = 0). 
(The equality is reached at £(T a ) = 0.) This means that if the third virial coefficient of the 
pion gas is negative at the temperatures T ^ 128 MeV, the model cannot describe the pion 
matter at these temperatures. Consideration of higher order terms of virial expansion would 
be necessary to obtain a satisfactory equation of state of the pion gas in this case. 

At T < T a the equation (|T4T ) has two solutions: n g and n x , n g < n a < n x (See Fig. ||). 
It is easy to check that the derivative of the pressure with respect to the density is negative 
at the point n x : 



dp(T, n) 



dn 



< 0. (19) 



n=n x 



This means that the root n x corresponds to an unstable phase, which cannot exist, because 
any small fluctuations of the density would destroy it. The solution n g corresponds to a 
stable phase which can be obviously interpreted as a pion gas. 

The interpretation of the second case is straightforward. The equation ( |I4"D has a unique 
solution at any temperature, which describes the gaseous phase of the pion matter. No PT 
takes place. 

The third case is the most interesting. The equation (|T4"D has solutions at any tem- 
perature. At low temperatures, T < T b , where T b satisfies the condition z(T b ,fi) = /(n&), 
the equation ([14]) has a single solution n g (n g < n a ), which can be obviously interpreted 
as the density of gaseous phase. As it is shown in Fig. |3], at intermediate temperatures 
T b < T < T a , T a is given by Eq . (|18D , two more solutions appear: n x {n a < n x < n b ) and n\ 
(ni > n b ). Similar to the second case, the root n x corresponds to an unstable phase, which 
cannot exist. In contrast, the points n g and n\ give a positive derivative of the pressure 
and correspond to (meta-)stable phases. The questions which of two phases is stable or 
metastable is resolved by the Gibbs criterion: the pressure of a stable phase is higher than 
that of metastable one. The point, when the two pressures are equal to each other 
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p(T pt , n g (T pt , /i)) = p(T pt , ni(T pt , //)) 



(20) 



is a PT point. It is reminiscent of the first order gas-liquid PT. At high temperatures T > T a 
the equation ( |T^| ) has an unique solution corresponding to this 'liquid' phase. 

In spite of its similarity to the gas-liquid PT, which is well known from molecular physics , 
there is an essential difference. In molecular physics, the stable liquid exists at temperatures 
below the PT point (T < T pt ), while high temperatures (T > T pt ) correspond to the gaseous 
phase. Here we observe the opposite situation: at low temperatures (T < T pt ) the pion 
matter only slightly differs from the ideal classical pion gas. The dense 'liquid' phase appears 
at high (T > T pt ) temperatures and thus can be called 'hot pion liquid'. 

There are two reasons that explain the qualitative difference in the behavior of the pion 
matter as compared with ordinary substances known from molecular physics: 

• In contrast to the number of molecules, the number of pions is not conserved and 
rapidly decreases as temperature falls down. So that the density of isotopically sym- 
metrical pion matter becomes so small at low temperatures that the interactions be- 
tween the particles plays almost no role, which corresponds to a perfect gas. 

• In molecular physics the dependence of the second virial coefficient on the temperature 
typically can be approximated by the formula [|J a>2(T) — (3 — where nonnegative 
constants a and (3 come, respectively, from attractive and repulsive parts of the in- 
termolecular potential. It is easily seen that at high temperatures the repulsive part 
becomes dominating and the virial coefficient becomes positive. In contrast, the sec- 
ond virial coefficient of the pion gas remains negative at high temperatures (see Fig.|l]), 
which indicates the attractive nature of 7T7r-interaction. 

As was mentioned above, we study the pion matter at chemical equilibrium. Thus the 
dense phase of pion matter is distinct from multipion droplets considered in Ref. [[LI]]. They 
may exist only at // > at low temperature and evaporate if the temperature exceeds some 
critical value, similarly to ordinary liquid. 

The dependence of the particle density on the temperature is shown in Fig. |j. The value 
of the free parameter £(T) = 0.21 was arbitrary chosen from the interval < £(T) < 1/3 
to provide an example. As is seen from the figure, the density of the gaseous phase is only 
slightly different from the ideal gas density indicating that virial expansion can be safely 
applied in this region. On the other hand, the density of the 'liquid' phase is about one 
order of magnitude higher than that of the pion gas, so one should not expect that the 
present model can provide quantitative description of the 'pion liquid'. It rather illustrates 
the features of PT at a qualitative level. 

At the same time, to estimate the upper bound of the PT temperature it is sufficient 
to analyze the properties of the gaseous phase. In fact, T pt cannot exceed the maximum 
temperature of the overheated gas T pt < T a (see Fig. [3]). The upper bound for T a is 
approached at £(T a ) — ► 1/3 and comprises 136 MeV. At £(T) > 1/3 the phase transition, as 
was mentioned above, disappears. That is, our model states that if a first order PT 
does exist in the isotopically symmetric equilibrium pion matter, it takes place 
at the temperature not higher than 136 MeV. 

A similar first order PT in hadron matter but at higher temperature (T = 190-^200 MeV) 



was predicted long ago within various models [12,13 1 . In Ref. [12] the hot and dense phase 



5 



was interpreted as quark-gluon plasma, i.e. the PT was identified with color deconfinement 
PT. This agrees with the estimations of the temperature of the deconfinement PT from 
lattice simulations of two-flavor QCD: T dec = 150 - 200 MeV 00 • 



In contrast, the present model predicts a PT at lower temperature, which cannot be 
identified with the deconfinement PT. This suggest the possibility that a new phase of hadron 
matter, distinct from the HG and QGP may exist at zero baryonic chemical potential in the 
temperature range 136 MeV ^$ T < Td ec - It is formed from colorless pions attracting to one 
another rather than from deconfined quarks and gluons and can be called 'hot pion liquid'. 

It is possible that the 'hot pion liquid' can be created in relativistic heavy ion collisions 
at RHIC and LHC energies, where the pions will play a dominant role at the final stage of 
the reaction. 



The authors are grateful to Larry McLerran for fruitful discussion and his interest to 
the paper and to Gephard Zeeb for reading the manuscript and usefull comments. The 
financial support of GSI and DAAD (Germany) is appreciated. The research described in 
this publication was made possible in part by Award No. UP1-2119 of the U.S. Civilian 
Research & Development Foundation for the Independent States of the Former Soviet Union 
(CRDF). 



6 



REFERENCES 



R. Dashen, S.K. Ma, H.J. Bernstein Phys. Rev. 187, 345 (1969). 
K. Huang, "Statistical mechanics" John Wiley & Sons (1963). 

J.E. Mayer and M. Goeppert Mayer, "Statistical mechanics" John Wiley & Sons (1977). 
W. Greiner, L. Neise and H. Stocker "Thermo dynamik und Statistische Mechanik" 
Verlag Harri Deutsch (1993). 
B. Hyams et aL, Nucl. Phys. B64, 134 (1973). 

S. Ishida, M. Ishida, H. Takahashi, T. Ishida, K. Takamatsu and T. Tsuru, particle," 
Prog. Theor. Phys. 95, 745 (1996) [hep-ph/9610325 . 



S. Ishida, T. Ishida, M. Ishida, K. Takamatsu and T. Tsuru, Prog. Theor. Phys. 98, 
1005 (1997). 

S. Ishida, |hep-ph/9905260 . 



A. Kostyuk, M. Gorenstein, H. Stocker and W. Greiner, |hep-ph/ 0004163 



A. Pais and G.E. Uhlenbeck, Phys. Rev. 116, 250 (1959). 

I. N. Mishustin and W. Greiner, J. Phys. G G19 (1993) L101. 

R. Hagedorn and J. Rafelski, CERN-TH-2947 Invited lecture given at Int. Symp. on 
Statistical Mechanics of Quarks and Hadrons, Bielefeld, Germany, Aug 24-31, 1980. 
J. Theis, G. Graebner, G. Buchwald, J. A. Maruhn, W. Greiner, H. Stocker and 
J. Polonyi, Phys. Rev. D28 (1983) 2286. 

T. Blum, L. Karkkainen, D. Toussaint and S. Gottlieb, Phys. Rev. D51 (1995) 5153 
||hep-lat/94100Tl . 



W. Greiner and D. H. Rischke, Phys. Rept. 264 (1996) 183. 



7 



FIGURES 
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FIG. 1. The dependence of the second virial coefficient of the ideal and real pion gas on the 
temperature. 
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FIG. 2. The dependence of the r.h.s. of the equation ( |l4|) f(n) on the variable n at // = 
for £(T) = and £(T) = 1/3 at T = T a = 128 MeV and T = T c = 136 MeV, respectively. The 
behavior of f(n) for £(T) < (^(T) > 1/3) is qualitatively the same as for ^(T) = (£(T) = 1/3). 
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FIG. 3. The dependence of the r.h.s. of the equation ( |l4| ) f(n) on the variable re at /j = for 
£(T) = 0.21 at the temperature T = T p t = 120 MeV. The behavior of f(n) for any £(T) from the 
interval < £(T) < 1/3 at any temperature is qualitatively the same. 
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FIG. 4. The dependence of the pion density on the temperature at chemical equilibrium \x = 
for £(T) = 0.21. The density of the ideal classical pion gas is shown for comparison. 
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